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Abstract 

> ' 

Q ■ We solve time-sliced path integrals of one-dimensional Coulomb system in an exact 

*y \ manner. In formulating path integrals, we make use of the Duru-Kleinert transformation 

with Fujikawa's gauge theoretical technique. Feynman kernels in the momentum represen- 
tation both for bound states and scattering states will be obtained with clear pole structure 
that explains the exactness of the path integral. The path integrals presented here can be, 
therefore, evaluated exactly by making use of Cauchy's integral theorem. 

1 Introduction 

<d : 

. One-dimensional system with the Coulomb potential, or the hydrogen atom in one space dimen- 

sion, may find its application not so much in practice!]] Nevertheless, after a work by LoudonpQ, 
. a number of publications on this apparently simple system have appeared [2] -[Hj to discuss the 

^ ' peculiarity of the singular potential. The strange result obtained by Loudon was the twofold 

degeneracy in the energy eigenvalues despite the one-dimensional nature of the system. Another 
interesting but quite controversial issue was the (non-)existence[3J |3J H] of the ground state of 
even parity with infinite binding energy. Exclusion of the even states with finite energy[5j and 
relativistic effect for the removal of such unphysical states 6 were also argued. Among these 
works, it seems that only a few attempts on analyses from the path integral formalism have so 
qq ' far been made. To our knowledge, only in Ref. [7] descriptions of functional-integral treatment, 

in addition to the functional-analytic approach, of the one-dimensional hydrogen atom can be 
found. The technique, based on the probabilistic measure, of the functional-integral employed 
in Ref. [7] is rather mathematical than that we wish to establish in this article. Furthermore the 
actual calculation of the functional-integral is worked out by making use of some identity that 
connects the Green function of the one-dimensional hydrogen atom to that of a four-dimensional 
harmonic oscillator. In other words, the path integral of the one-dimensional hydrogen atom 
is not still solved explicitly. It cannot be, therefore, sufficient for our purpose to understand 
the exactness of the path integral for one-dimensional hydrogen atom by such calculation. As 
stated in Ref. [7J , there may exist some natural prescription for the boundary condition of wave 
functions at the origin required from the consistency of the path integral. We should be able 
to find such structure by constructing path integrals in the time-sliced form keeping connection 
with the operator formalism. Therefore, from theoretical as well as pedagogical points of view, 
it will be useful to develop another path integral technique for solving such a singular system. 

Path integral of the hydrogen atom was first successfully formulated by Duru and Kleinert[Sl 
9 by making use of the so called Duru-Kleinert (DK) transformation. The DK transformation 
consists of two fundamental ingredients: reparametrization of the time in the path integral and 
making use of the "square root of the coordinates" by means of a Kustaanheimo-Stiefel(KS) 
transformation [TUj . Many attempts, mainly on the use of KS transformation, by which we can 
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1 After the completion of this work, the author has been made aware of experimental studies on the exciton 
spectra 24 under the high magnetic field. It is shown in Ref. 25 that an exciton under extremely high magnetic 
field can be well approximated by the one-dimensional Coulomb system. 



1 



write the Coulomb path integral in terms of oscillator coordinates, with or without time-slicing, 
have been made to clarify the approach by DK transformation[Il]-[22]. It was Fujikawa [23] 
who first observed that the essence of the DK transformation is the special choice of the gauge 
fixing condition for a system with invariance under reparametrization of the time. He proposed 
a new approach to formulate a Coulomb path integral in an exactly solvable way from the view 
point of the Jacobi's principle of least action. He also noticed that the use of KS transformation 
in the Coulomb path integral is rather technical and demonstrated the exact solubility of the 
Coulomb path integral in three dimension in terms of the parabolic coordinates. 

The hydrogen atom in one space dimension seems to be rather exceptional in the path 
integral formalism. Although many examples with detailed proofs are found in the textbook[22j 
by Kleinert, the space dimension Dc for the path integral of Coulomb potential in chapter 14 
of the Ref. 22J is assumed to be Dc > 2. In particular, if we set Dc = 1 in Eq. (14.93) of 
Ref. [22], we obtain Do = as the space dimension of the corresponding, or DK-equivalent, 
oscillator coordinates. In fact, there does not seem to be known any KS-type transformation 
for one-dimensional Coulomb system excepting the one utilized in Ref. 7] as mentioned above. 
Motivated by these observations, we adopt in this article Fujikawa's formulation to find path 
integral descriptions for one-dimensional Coulomb potential V(x) = —a/r, where r = \[x^ . Use 
of the KS transformation will be thus avoided in our formulation. It was the problem of operator 
ordering, arising from nonlinear canonical transformations such as KS transformation, in path 
integrals that forced some authors to criticize the DK method with KS transformation(see 
Ref. [20] . [2"2"| and references therein). Our formulation of time-sliced path integrals in this 
article will be built basing on a symmetric Hamiltonian operator that has definite rule of 
operator ordering, thereby being free from the ordering problem. The effect of time-slicing 
from the formal continuum path integral will also be problematic if we rely on such formal 
derivation of path integrals. In our formalism, the basic ingredient is the short-time kernel 
that describes an infinitesimal time evolution. Feynman kernels of finite time evolution will 
be constructed by multiplication of such infinitesimal kernels. Without this process it will be 
impossible to understand the mechanism of exactness in path integrals considered here. 

We organize this paper as follows: in section 2 we will briefly explain Fujikawa's method then 
formulate immediately a path integral with a negative fixed energy in the Euclidean formalism. 
For the case of positive fixed energy, the Euclidean technique does not work well. Hence it 
will be treated separately for the usual Feynman kernel by setting the time parameter to be 
real. For both of these path integrals we shall develop a novel procedure in multiplication of 
infinitesimal kernels in the momentum representation. The use of the momentum space wave 
functions is due to the form of the symmetric Hamiltonian we set up. This will become clear 
soon in the next section. Section 3 will be devoted to the construction of the physical states 
from the solution of path integrals in section 2. In Fujikawa's procedure this is the important 
step for solving the problem in terms of path integral. This will be achieved by analyses of the 
resolvent operator. Conclusions and remarks will be found in the final section. 

2 Path Integrals of the Free Particle with a Fixed Energy 

For a system, whose Lagrangian being given by L or i g = x 2 /2 — V{x), we consider another 
Lagrangian with reparametrization invariance, given by 



for a fixed energy E, according to the procedure introduced by Fujikawa in Ref. [23] . The 
invariance of this new system under a time reparametrization requires a gauge fixing upon 
quantization. As was shown by Fujikawa, by choosing an appropriate gauge condition, we can 
formulate a path integral that is equivalent to the DK transform of the path integral for the 
original system. 

The Lagrangian ([T]) yields the canonical momentum, given by 




(1) 
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that results in a constraint 

0o := \p 2 + V{x) -E = 0. (3) 
In addition, due to the fact H = px — L = 0, we need to test Dirac's total Hamiltonian, given 

by 

H T = u(x,p)<t> ttO, (4) 

in which the Lagrange multiplier u(x,p) cannot be determined from the equation of motion at 
all. Choosing this rather arbitrary degree by hand exactly corresponds to the choice of gauge 
fixing condition. The convenient way for us is setting u(x,p) — r to get 

Hi = 2 V P 2 ~ a ~~ (5) 

for the one-dimensional Coulomb potential. However, there is no restriction in application of this 
procedure and gauge fixing to systems other than the Coulomb system. We may, for example, 
enjoy this method even for the free particle by setting a = in the above. Furthermore, if we 
solve path integrals of the free particle in this formalism, the solution of corresponding ones for 
the Coulomb system is immediate, since the difference in the Hamiltonian Hi of a ^ from 
that of a — is merely a constant. We therefore first examine the quantization and its path 
integral of the total Hamiltonian ([5]) for the case of a = in the following. 



2.1 Path integral for the negative energy 



Quantization of the system can be achieved by adopting a unitary representation of the canonical 
commutation relation(CCR) 

[x,p]=i (6) 

with the constraint Q as a physical state condition. Since the Hamiltonian ifp involves a term 
that consists of r multiplied by p 2 , we need to specify an operator ordering upon quantization. 
For brevity we choose a symmetric product for this term to consider the operator Hamiltonian 



1 



-prp 



Ef. 



(7) 



To elucidate the property of this Hamiltonian, we here put E — —uj 2 /2 and evaluate a Euclidean 
kernel(a Feynman kernel with imaginary time) in the momentum representation: 



K{pb,n\p a ,T a ) = (p b \e 



-Hs(Tb—T a 



\Pa), 



(8) 



where \p) designates the momentum eigenvector. Since the Hamiltonian ([7]) is linear in f, the 
integration over x in a path integral can be easily performed. This is the reason why we take 
the momentum(p) representation instead of the coordinate^) representation. 

The first step for the evaluation of the kernel is the calculation of an infinitesimal kernel 



(p\ (l-eHs) \p') 



dx 

^ 6XP 



-i(p — p')x 



j(pp 



(9) 



where e = (r& — r a )/N > is taken to be uje -c 1. We split the integration domain at x = 
and put x —x for negative x to find 



(p\ (l - eH s ) \p>) = exp [- {i(p 



' ~ P') a + 2 (pp' 



(10) 



To ensure the convergence of this integration, a condition, pp' + lo 2 > 0, is necessary. Although 
it is not explicit, this condition is equivalent to p 2 + lo 2 > 0(p' 2 + uo 2 > 0) that holds always 
for real p(p'). To check this we rewrite pp' as p 2 — p(p — p') in the above to make a change of 
variable from x to (1 + iepa/2)x. By taking the Jacobian into account, we find that the right 
hand side of Eq. pH|) is equivalent to 



E 

<T=±1 



'dx 
2^ 



exp 



-pa 



{i(p - p')a + | (p 2 + lj 2 )}. 



(11) 
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up to 0(e) in the exponent. Note that the additional term from the Jacobian precisely matches 
to the one we need when rewriting the symmetric product pfp as p 2 f + [f,p\. The condition 
for convergence of the integration now reads p 2 + uj 2 > for this new expression. This proves 
the above statement. We can therefore proceed to carry out the integration with respect to x 
to find 



(p| (l - eiTs) \p) 



1 

2nd 



1 



p — p' — ie (pp' + uj 2 ) /2 p — p' + ie (pp' + uj 2 ) /2 



(12) 



Once we obtain the expression (|12p , we can regard it as a function of two complex variables p 
and p' . As a function of p, each term in the right hand side has a simple pole at 



P : 



p' + ieuj 2 /2 p' - ieuj 2 /2 
l-iep'/2 ' l + iep'/2 ' 



(13) 



respectively. As a confirmation, we observe that the real part of pp' + ui 2 is positive for both 
of these poles if we set p' to be real in this evaluation. It will be interesting to see that the 
situation is quite similar to the one for an expression of a delta function 



8(p-p) 



1 



1 



1 



2m \p — p' — i0 + p — p' + i0+ 



(14) 



where 0+ designates a positive infinitesimal. We may often utilize this notation hereafter. 

The second step for us to proceed is the multiplication of the infinitesimal kernels. But 
before doing this we observe 



p - p' T i\ (pp' + ^ 2 ) = (P - P') cosh(we/2) =F - (pp' + u 2 ) sinh(we/2) + 0(e 2 ) 
2 uj 

so that we can imagine the kernel for a finite imaginary time r(r = t/, — r a ) to be 



K(pb,n\p a ,T a ) 



E - 



<T=±1 



2?ri 



(Pb - Pa) cosh(cjr/2) - i— (p b p a + uj 2 ) sinh(wr/2) 



(15) 



(16) 



Let us prove this by explicitly evaluating a product of infinitesimal kernels. To this aim we 
recall the analogy to the delta function above again. The first term of the infinitesimal kernel 
(fT2"|) possesses a simple pole of positive (negative) imaginary part as a function of p(p') while the 
situation switches to the opposite in the second term. We can therefore make use of the Cauchy's 
integral formula, just like the case we will do in the multiplication of two delta functions, in 
carrying out the integration 

(pa| (l - eH s y |po) = / d Pl (p 2 \ (l - eH s ) \ Pl )( Pl \ (l - e£ s ) |p ), (17) 



to obtain 



(pa| (l-ei?s) 2 |Po) = ^- 



This suggests us to suppose 



(P2 - Po){l + (uJe/2) 2 } - iuje— (p 2 p Q + uj 2 ) 



(P2 - Po){l + (we/2) 2 } + iuje- (p 2 p + uj 2 ) 



(18) 



( Pj \ (1-eHs) bo) = 



2iri 



1 1 

(Pj - Po)C 3 Sj (pjpo + uj 2 ) (pj - p )Cj + -Sj (pjp + uj 2 ) 

UJ UJ 



(19) 
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with initial conditions, Cq — 1 and 5*0 = 0+. Then, by repeating the same procedure above, we 
find that one more multiplication by a short-time kernel (pj+i\ ^1 — eHsj \pj) to this expression 
results in 



(p i+1 | (I - eH s ) 1+1 \ PQ ) = ^- 



x 



E 



2wi 

a (20) 



<t i:i (Pj+i -Po)(Cj +ueSj/2) - i — (veCj/2 + Sj) (Pj+iPo + ^ 2 ) 



By comparing Eq. ([2"U|) with Eq. (fm|) . we obtain a set of recurrence formulas 

Cj+i = Cj + -^Sj, Sj+i = — Cj + Sj. (21) 

The solution of the above can be found immediately to be 

Cj = cosh(jwe/2), Sj = sinh(jwe/2) (22) 



to convince us that our conjecture in Eq. (|16|) holds. 

Eigenvalues and eigenfunctions of the Hamiltonian H§ can be extracted from the kernel in 
Eq. (rTB|) . To see this we rewrite the first term of K(pb, Tb\p a , i~ a ) as follows 

1 1 



2?Fl (Pb - Pa) cosh(wr/2) - — (pbPa + u 2 ) sinh(wr/2) 

U) 



1 



(23) 



7r e" T / 2 (cj + ipb){uJ - ipa) - e UT / 2 (u - ip b ){uj + ip a ) 

J± l _ \p e -(n+l/2) WT ( u-iPb Y ( u + iPa 

7T (ui + ipb)(ui - ip a ) £^ Q \uj + ipb) \uJ-ipa 

In the same way, the second term is also expanded in a series as 
1 1 

27T2 4 

iPb ~ Pa) cosh(wr/2) + - (pbPa + UJ 2 ) sinh(wr/2) 

UJ 

J± 1 e -(n+l/2) WT f u + ipb X' 1 ( u-ip a 

7r (u) - ip b )(uJ + ip a ) \UJ-ipbJ \uj + ip a/ 

We thus obtain another expression of the kernel in the series of eigenfunctions, given by 

oo 

K(pb,T b \p a ,r a ) = £ e -(«+i/2)- {^(p b )^+Hp a r + ^-Hp b )^-Hp a r}, (25) 

?i=0 



(24) 



where ijj^ip) for n = 0, 1, 2, ... , defined as 



^ ±) W = V-4-(^J . (26) 

V 7TW±Zp \LU ±ip / 

are eigenfunctions belonging to a common eigenvalue An (E) = (n + l/2)u in the momentum 
representation. From these manipulation we observe the twofold degeneracy in all eigenvalues 
of the Hamiltonian. Eigenfunctions in the position representation can be obtained from ^„ (p) 
as 

^\x)= [°° ^Le^^ip). (27) 



— DO 



Since i/>n + ' ) (p) (V>n ^(p)) has a pole of order n + 1 at p = iw(p = — iw), the eigenfunction 
(a;) {iph \x)) has support only on positive(negative) x. Hence we find 

(_l\n jn 

i> { n +) (x) = ( -^9(x)V2^ — [e^to + iu)*] 

nl d P P =iu> (28) 

= (-l) n V2^9(x)e- UJX L n (2ujx), 



where L n (z) designates a Laguerre polynomial of the n-th order and 0(x) is the step function. 
In the same manner, by calculating the residue at p — —ilu, we obtain 



^ _) 0) = (-l) n V2^6(-x)e u ' x L n (-2ujx). 

In short, we may write 

^ { n\ x ) = (-l) n V2^9(±x)e- ur L n (2Lur). 
In this way, by writing the kernel (xb\e~ HsT \x a ) as K(xb, Tb\x a , r a ), we obtain 



(29) 
(30) 



OO 

K(x b ,n\x a ,r a ) = Y^z~ (n+1/2)ulT {^\x h )^+\x a )*+^-\x b )^-\x a y), (31) 



n=0 



corresponding to Eq. 

We have thus observed that a time-sliced path integral, 



f N dp N 

J n 2^ n dx ? ex p 



4=0 3=1 



JV 



k=0 



N 



i^2pk(xk+i - Xk) - ^ (pkPk-i + lu 2 ) r k 



fc=i 



(32) 



that describes a kernel K(x b , T b \x a , r ) in the limit iV — ► oo, can be evaluated exactly by making 
multiple use of Cauchy's integral theorem. Therefore we can regard the model considered here 
an interesting example of an exact path integral that can be carried out without use of the 
Gaussian identity. 



2.2 Path integral for the positive energy 

So far our consideration has been restricted to the Hamiltonian Hg = pfp/2 — Ef with a negative 
energy E = —uj 2 /2. It will be also interesting to see how the path integral of this system with 
a positive energy E = k 2 /2(k > 0) can be solved. Hence our next target to be examined here 
is the time-sliced path integral 



N 



n n dx i ex p 



N 



2tt 

i=0 3 = 1 



>2_ /',(.'•,. i .r ; j 

3=0 



. N 
It x ^ 

2"^ 

3 = 1 



PjPj-l - k2 



(33) 



Unfortunately, however, we cannot define its Euclidean version without restricting the range 
of values for p as well as for the imaginary time r — —it due to the lack of positivity of the 
Hamiltonian in the term proportional to p 2 — k 2 . We shall therefore briefly point out some 
properties of the Euclidean path integral for this system at the end of this subsection. 

We begin again with evaluation of the kernel in the momentum representation. By adding a 
positive infinitesimal to the imaginary part of k 2 , we can perform integration over x in Eq. (|33[) 
to obtain a short-time kernel 



(Pl( 



1 - ieH s 



1 



1 



(P -p') + e (pp' - k 2 ) j2 {p ~p')~e (pp 1 - k 2 ) j2 



(34) 



Thanks to the regularization factor, poles of each term in this kernel are shifted off to the 
imaginary axis from the real line. For instance, as a function of p, the kernel (p\ ^1 — ieHs^j \p') 



has poles at 



P 



p' + ek 2 /2 



and 



p' - ek 2 /2 



(35) 



1 + ep'/2 " 1 - ep'/2 ' 

whose imaginary part being positive and negative, respectively. Once we recognize this pole 
structure of the infinitesimal kernel, we can immediately carry out integrations over p for 
multiplication of short-time kernels by making use of the similar technique we have employed 
for the case of negative energy. We thus omit showing here the detail of the procedure for 
finding a Feynman kernel 



K(pb,tb\p a ,t a ) 



1 



2wi 

<7=±1 



(Pb - Pa) cosh(fci/2) + - (p b p a - k 2 ) sinh(fci/2) 
k 



(36) 



G 



in the momentum representation for a finite time t = % — t a > 0. Note that the regulator for a 
short-time kernel is included even in this expression as infinitesimal shifts T*0+ for a = ±1 in 
the denominator. 

If we set p = kcoth.8 or p — fctanh# corresponding to the sign of p 2 — k 2 , it becomes easy 
to decompose the above kernel into an integration over the continuous eigenvalue of Hg. At 
the outset, we set p a — fctanh# a and p = fctanh^. Then we rewrite each term in the kernel 
to obtain 



1 cosh Oh cosh 9„ 



( Pb - Pa ) cosh(fci/2) ± 1 ( PbPa - k 2 ) sinh(fci/2) k Sinh ^ b 6a T kt l^ 
A useful formula for our aim here is an integral representation 

1 r°° 1 



(37) 



sinh x 



= 2i / dz- ^e 2lxz , (38) 

1 + e 2 ^ v ; 



in which a condition, — 7r < Q(x) < 0, for the convergence being necessary. We here remember 
infinitesimal shifts =F*0+ in the Feynman kernel (|36[) to regard ±(9 — 9 a ) — kt/2 as x in Eq. 



Then by noticing that \Jk 2 — p 2 = kj cosh 9 for p — k tanh 9 , we can express the kernel Eq. (|36|) 

as 

k 1 

K(p b ,t b \p a ,t a ) = - 

77 V( k ~Pb)( k -Pa) 

(39) 



X V / due- iHv 



1 





+ p 6 \ , 


(k 


+ Pa\ 


(r 


-PbJ 


\k 


~Pa J 



^ / 1 

for a combination of p a and p b that satisfy [p | < k and < A:. 

In the same way, for p a and p b satisfying \p a \ > k, \p b \ > k as well as p a p b — k 2 > 0, we may 
set p a = k coth # a and p b = k coth to find 

1 sinh 8 b sinh 9 a 



(p 6 - p a ) cosh(fc</2) ± i(p bPa - sinh(fci/2) fcsinh (^ *« T fc< /2) 

in which 9 b — 8 a should be regarded as 9 b — 9 a ± i0 + corresponding to the sign in front of kt/2. 
We then utilize the formula (f3"8"|) again to obtain 

k 1 

K(p b ,t b \p a ,t a ) = 



* V(Pi-k 2 )(p 2 a-k 2 ) 
^— ' / 1 + r 27 



X 

<J=±1 



Ph-KJ \Pa 



(41) 



for this case. Here we have made a change of variable from v to — v and rewrite 1/(1 + e~ 27UJ ) 
as e 2TW /(l + e 27r ") to reach this expression. Although there exist other possible combinations 
of values for p a and pt>, it will be enough to observe Eq. (j3"9"|) as well as Eq. flTj) for finding 
eigenvalues and eigenfunctions, ^ (p) = (p|A;±), of H$ in the momentum representation. 
Explicitly, wave functions are given by 



~ (±) 1 (k +p)(-!± 2 ^)/ 2 

^ b)= ^(1 + e^) (jfc_ p )(^)/2 (42) 



so that | A; ±) to be eigenvectors of Hs with an eigenvalue A = kv\ 

H s \X;±) = ku\X;±). (43) 



Note that the regularization for the short-time kernel Eq. (l34[) can now be interpreted as the 
prescription for avoiding branch points on the complex p-plane. In the region \p\ < k, p is 
shifted with an amount of =r0+ in 4>^{p) by the effect of the regulator. On the other hand, 
when p is in regions where \p\ > k, the shifts occur in the opposite way. This specifies the 
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way to go around the branch points at p = ±fc. When p moves along the real axis, a detour 
that goes lower and upper side of p — ±k is made for ^p^\p) and for -0^ \p), respectively. 
Due to these singularities of branch points, discontinuities of wave functions as functions of a 
real variable p at p — ±fc are unavoidable; if we consider the analytical continuation of the 
eigenfunction (p) in the region \p\ < k to that of p > k, for instance, the wave function gets 
a factor e - l7T / 2 + 7r ' y multiplied in addition to a switch from k — p to p — k in the denominator. 
Nevertheless, we can clarify the completeness of the eigenf unctions simply by setting t = in 
each expression, Eq. (|39[) . (|4ip . and in other similar ones, of the Feynman kernel. If we keep in 
mind the discontinuity of eigenfunctions at branch points, we can also easily check the following 
normalization and the orthogonality relations of |A; ±), 

(\;a\\';a') = 5 a , a ,6(\-\'). (44) 

We have thus completed the calculation of the Feynman kernel in the momentum representation. 

Turning now to the original time-sliced path integral of Eq. (|33p that defines a Feynman 
kernel in the position diagonal representation, we make use of the integral representation of the 
confluent hypergeometric, or Kummer, function to obtain 



(k+p) 



(-l+2ii/)/2 



(fc_ p )(l+2 



iu)/2 poo -i 

— = J dx e-^'^Fi- - iv, 1, -2ikx), 



(45) 



where the shift of p by an amount of — z0+ must be understood. With the help of this formula, 
we can find the eigenfunction, (x\\;+) — ip^\x), for an eigenvalue A = kv in the position 
representation as 

^\x) = 6(x)J——e* kx F(- - iv, 1, -2ikx). (46) 
V cosh7w 2 

In the same way, (p|A; — ) can be transformed to {x\\; — ) = ip^ \ x ) as 



^-\x) = 9(-x)J -^—e- ikx F{\ - iv, 1, 2ikx). (47) 

" COSh7Ti/ 2 



These results enable us to transform the Feynman kernel K(pb, tb\p a , t a ) in the momentum 
representation to the one in the position representation. Writing the Fourier transform of 
K(pb, tb\p a , t a ) as K(xb,h\x a , t a ),we finally observe that the time-sliced path integral of Eg. ([55)) 
can be carried out exactly to yield 

/oo 
dX e- lXt { ^ +) (x b )ip{ +) (x a )*+ {x b )^x ] (»«)* } ■ (48) 
-oo 

It will be obvious that the essence of the exact calculation in the time-sliced path integral is 
again the pole structure of the short-time kernel Eq. (|34|) similar to the one we have found in 
Eq. (fT2")) for the case of a negative energy E — —uj 2 /2. 

As we have stated above at the beginning of this subsection, the Euclidean path integral of 
the Hamiltonian H§ with a positive energy E = k 2 /2 cannot be freely defined. Nevertheless, we 
may try to set t — — ir (Wick-rotation) in expressions, Eq. ([39| and (|4"Tj) . of the Feynman kernel 
K (pb, t \p a , t a ) in the momentum representation. If we make this change in Eq. (|41[) , we find 
a condition, < kr < 2tt, for the convergence of the integration over v. On the other hand, 
Eq. (f3"5|) can be Wick-rotated only for the case —2tt < fcr < holds. Hence a kernel, though 
existed, for the imaginary time r cannot be extended its domain of definition into other regions 
of momentum variables. For this reason we are restricted to make considerations on Feynman 
kernels with t in the real domain. 



3 Physical States 

We have clarified in the previous section that path integrals of the Hamiltonian Hg for Feynman 
kernels with both positive and negative energies can be performed in an exact manner. The 
Feynman kernel of Hs is, however, the one for the unphysical system in a specific gauge. To get 
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information of the original system, we need to take the physical state condition ([!]) into account. 
This can be achieved by considering the resolvent operator (H — E) , where H expresses the 
Hamiltonian of the original system without gauge invariance. In this section, we first examine 
the resolvent operator for a free particle with a negative energy E = —uj 2 /2, then generalize 
the method to fit the Coulomb system by making use of the result in the previous section. 
Establishing the prescription to find the resolvent operator from the Feynman kernel in these 
examples, we eventually solve the scattering states of the one-dimensional Coulomb system. 

For the case of a free particle governed by the Hamiltonian H = p 2 /2, we first define 
non-symmetric operator = (H — E)r(E = —lo 2 /2) to observe 



(h-E) 1 =fH R \ (49) 

To evaluate the right hand side of the above in terms of the path integral, we here invent 
a useful relation that connects H-r. to the symmetric one Ha,. To this aim we introduce an 
operator a by a — x/f, which can be defined excepting the origin x = 0, to express the sign of 
the position operator. Note that it is also possible to introduce a by the commutator [f,p] = ia 
as we have already done implicitly in obtaining Eq. (|ll[) . As an important fact, it should be 
remarked here that a commutes with Hg = pfp/2 — Er, i.e. aHg = Hga. It will be clear that 
the operator a is responsible for the supplementary quantum number a in the eigenfunctions of 
the Hamiltonian Hg. Needless to say, it is not a good quantum number for Hamiltonians other 
than Hs- Nevertheless it will be useful for finding eigenfunctions of the Coulomb Hamiltonian 
by restricting the domain of x to one of a half line corresponding to the value of a. Keeping 
these in mind, we now define another operator u> ~ ip and its Hermitian conjugate to find 

H R = (ujTip)(H s ±^(LUTwr 1 ■ (50) 

Though both of these formulas are equally useful for our purpose, we here choose the one of 
lu — ip for the formulation below. 

By writing V'i ± ^(p)(n = 0, 1, 2, . . . ) in Eq. ([2"6"|) as (p\n; ±), we denote right eigenvectors of 
Hs as \n; ±). Then we define IV'r"'^) = (w — ip)\n; ±) to find 

e~^>R ;+) > = e-(" +1 ^ T |4" ;+) ), e~"« T \^- } ) = e'™^ ; " } ). (51) 

: eigei 

gates, defined by ('0 R n ' ± ' l | = (n;±|(w — ip) -1 , to yield 



Completeness of the eigenvectors of H$ will be expressed in terms of IV'r ) anc ^ their conju 



Sandwiching between (pf,\ and \p a ) this expression after multiplying f from the left, we obtain 
(p b \fe- 6 ^\p a ) 

CO J' C Q^j 

Evaluation of the operator f in the above is done as follows: we may write f = xa in which x 
has clear definition in the momentum representation, given by (p\x = id(p\/dp, our task reduces 
to calculate (p\d'\ip R ' crS> ) for a = ±1. By noticing [<r,p] = 2iS(x), where the delta function of x 
being defined by 

1 f°° 

m = 7T d P' eWX ( 54 ) 



2tt 



— DO 



so that 



(pI*(*) = ^: / <ii>'(i>'\- O'V pVI- i->- r >) 



we rewrite ip\^\^^^) as 



(p\ {(w - ip)a- i[a,p]} \n;a) = a{u - ip)(p\n;<r) + - / dp' (p'\n;o). (56) 

7T 



oo 

' /—'I 
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Since the second term above is independent of p to be destroyed by differentiation with respect 
to p, only the first term contributes to (p\f\ip^' a ' > } to yield 

(p\r\^R ±] ) = ±^ {(w - ip)(p\n; ±)} . (57) 

Remembering the definition of wave functions ipn (p) = (p\n;±) in Eq. ()26|) . we may write 
explicitly as 

. , 

7T L0 Z + p z \UJ — ip / 

Here it must be remarked that (pl^l ' ) vanishes for n = thereby disappearing from the 
series in Eq. (|53p . Combining the result above with those of (^r 1 \Pa), we find 

( Pb \fe-^\p a ) =f]e-"" T - -, „ 2 ™„ 2 , 

5i k{u 2 +pD{^+pd (5g) 

U) — ipb U> + ip a \" ( W + iPb W — iPa * 



7T LO^+p^ \UJ + ip y 

( n ._v Ilo 2nu! ( lo + ip^ " 



^LO + iph LO — ip a J \U) — ipb LO + ip a / 

We are now at the point to observe, by integrating the above with respect to r, that 
.1 2 to 2 

(P&l "7v ~bo> = 7-0-; 2\ S (Pb ~Pa) -f-~2~i 2\( 2 i 2~T ( 60 ) 

H-E (u 2 +p£) ir {lo 2 + P 2 )(lo 2 + p 2 a ) 

for E = —lo 2 /2. It is precisely the first term of the right hand side that should be defined by the 
left hand side. Therefore the second term above indicates a contradiction in our manipulation. 
The origin of this obvious disagreement is the disappearance of (p^IV^r ) from the sum. A 
quick recovery for this discrepancy will be possible if we keep (p|f ) until we carry out 

the integration over r, by replacing it with e~ u ^ T {p\r\^->^' ), thereupon letting v — ► 0. To 
be consistent, however, we should examine this phenomenon in detail otherwise we will not 
be successful in establishing the method to define the operator fH^ . Though interesting, it 
will not be so useful to look more closely at the resolvent operator with a negative energy for 
the free particle. Hence we terminate the analysis on this system here by learning that the 
definition of the resolvent operator by its action on each eigenfunction of the non-symmetric 
operator fails due to the disappearance of an eigenfunction from the complete set. 

We now switch the Hamiltonian to that of the Coulomb system by changing to Hr — a 
in the above argument. Observation of the failure in the naive approach above by integration 
over r term by term on each eigenfunction of H-r suggests us to try dealing with each term in 
the Euclidean kernel 

(p6|e-<* 8 - a > T |Pa) 

= JL Z. ( 61 ) 

a=±i 2lTl {Pb - Pa) cosh(wr/2) - i— (p b p a + lo 2 ) sinh(wr/2) 

LO 

as a whole. To this aim, let us set pb — lo tan 6b and p a = lo tan a then rewrite the above in 
terms of these new variables. Firstly, we write 



(Pb - Pa) cosh(^T/2) =F — (pbPa + lo 2 ) sinh(wr/2) 

(62) 

% -9 a T W2) 



to find 



cos 9b cos 6 a 
( Pb \re-^ a ^\p a ) 



COS 2 b COS 2 6 a d e -i(8b-8a-i™T/2)+ a T 

<J = ±1 



2ttw 2 d0 b sm{0 b -6 a - icrLor/2) (63) 
cos 2 Ob cos 2 a 



2ttlo 2 



sm 2 (0 b -9a- ilot/2) ' sm. 2 {0 b - a + ilot/2) 
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Since we are dealing with the case a > 0(attractive Coulomb potential), both of the terms in 
the above monotonically tends to if we let r — > — oo. Although the path integral for the 
Euclidean kernel has been performed only for r > 0, the result can be regarded as a function of 
a complex variable r, besides being a function of pb and p a . The above mentioned behavior in 
the asymptotic region, r — ► — oo, naturally suggests us to integrate Eq. (|63p with respect to r 
form to — oo for a definition of f(HR — a) -1 . However, when the interaction is turned off by 
setting a = 0, the behavior of the first term becomes wrong. This will be confirmed easily by 
looking at the second line in the above for a = +1 and t —> — oo. Such a discrepancy must be 
avoided for the consistency of the formalism. We therefore choose the domain for integration 
over r from to +oo for a = +1 and from to — oo for a = — 1 with the assumption lu > a for 
a while. Note that, by this new prescription, we can explain the inconsistency found in Eq. ()60p 
above. The undesirable term is the consequence of the non-vanishing, though finite, behavior 
of the second term in Eq. (|63p when integrated over r from to +oo(o! = should be set for 
the free particle). 

We here utilize the integral representation for a hypergeometric function to carry out the 
integration of Eq. (j6"3"|) over r, thereby obtaining 

/ d JL =_if ^f(2,lTa,2Ta;e-*- 9 -'), (64) 

Jo sin (#6 — 6 a — iujT/2) lu 1 =p a 

where d = a/a; and F(a,b, c; z) denotes a hypergeometric function. Here we should recall the 
initial condition, So = 0+, imposed in the derivation of the Euclidean kernel. We have chosen 
0+ instead of to make the kernel reproduce the delta function 5 (pi, ~p a ) even after taking the 
limit t — ► 0. This prescription now takes effect in the above by multiplying e T °+ to e - 2l (0b~8 a ) 
depending on the sign in front of a. In view of this prescription, we should rewrite one of the 
above as 

-F(2,l + d,2 + d;z) = — ^r-F(2,l-a,2 -&;-) + -?^r(-z)- & (65) 



1 + a 1 — a z 

for z = e -M e b-0a+iO + ) t De useful as a hypergeometric series. Taking these into account, we 
obtain 

' / -W(«»-9.-<0 + ) F ( 2j 1 _ d; 2 „ & . e -2i(6 b -9 a -i0 + )j 



1 -a 



(66) 

+e 2i(9 b -9 a +io + ) F ( 2> 1 _ fij 2 _ & . e Me b -e a+l o + )^ 



simva 



_ e -2i(e b -e a +io+) 



for the Hamiltonian Hq = p 2 /2 — a/f with E — —lo 2 /2. Here use has been made of the fact 
that the effect of our prescription of adding ±0+ in the denominator of the kernel is equivalent 
to regarding E = —uj 2 /2 as to be associated with +10+. Series expansions of hypergeometric 
functions in the above will then bring us 



(Pb I A ^ — \Pa 



H c -E-i0 + a n (u 2 + pI)(lu 2 + p 2 a ) 

n ( / LP - ipi, LU + ip a \ " / id ~j~ ipi, CO ipg 
^ n - a \ \u> + ipi, u> - ip a J \ lu - ip b 10 + ip a 



(67) 



2nia ( lu — ipi, lu + ip a 



— a 



1 — e 27rm y u) + ipb lu — ip a 

Clearly, there exist poles of the Green function (|67|) at E — E n 

lu 2 a 



E n = —£, u> n = -, (68) 



n 
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for n = 1,2, 3,..., corresponding to bound states with twofold degeneracy for each energy 
eigenvalue. By approximating to(E) 
n-th eigenvalue, we obtain 



V-2E as uj{E) ~ uj n + (E n 
'lu =fip b Lu± ipa 



E)/u> n in the vicinity of 



7T (lu 2 + pD(u 2 + pi) n-a 
1 2 /a\3 1 



E n — E 



7r V nJ 



LU±ip b LUT ip a 

LU n T »Pb ± »Pa 



(69) 



to extract the pole parts of the Green function as 

i , . A i 



(P6 



H c - E - iO. 



-\Pa) 



n— 1 



+ ^(p 6 )$^(p o r} (70) 



where 



$£(p) 



7r (n) 



3/2 



1 



+ P \^n± IP ! 

We have thus clarified the pole structure of the Green function. Apart from the continuous 
spectrum, poles of the Green function precisely describe the zeros of He — E. The gauge 
fixing condition, given by Ht ~ 0, in our formulation of the corresponding classical system 
is proportional to He — E with an almost non-vanishing factor r. Therefore we may regard 
enforcing He — E = in quantum mechanics as the implementation of the gauge fixing condition 
that picks up physical or admissible states from the Hilbert space of the system with a gauge 
degree. 

Another Green function that reflects different boundary condition is obtained by changing 
the prescription from E — > E + i0 + to E — > E — i0+, instead. For this case, poles appear at u> = 
—u> n . These cannot be regarded as bound states since > is needed in setting E = —oj 2 /2. 
Wave functions for continuous eigenvalues can be found by subtracting (He — E + i0+) _1 from 
(He — E — i0 + )^ 1 . It will be useful here to consider in the region E = k 2 /2 > to find 
the continuous eigenfunctions. The transition from a negative energy to a positive one can be 
most easily achieved by setting lu — > — ik, 9b, a — * *#&,a as well as writing a/k as a again. The 
expression for E = —lu 2 /2 in Eq. (f6"3"|) then becomes 

(pblre-^-^lpa) 



cosh 2 9b cosh 2 9 a 
2rf 



(72) 



_simV(0 b - 6 a - ikt/2) sinli (6> 6 - 9 a + ikt/2)_ 
This expression will be integrated over t from to +oo(— oo) for the first(second) term to yield 
1 k 2 




IT (k 2 - P 2 )(k 2 

k + Pb k - p a 

Pbk+p a 



Pi) 



k - p b k + p a 
k + Pb k - p a 



(73) 



in which a = a/k has a negative imaginary part due to the boundary condition of the Green 
function: k — > k + i0 + for E — > E + i0 + . Here we have made use of the formula Eq . ([38| again in 
obtaining the above. If we change the boundary condition by replacing this with E — * E — i0+, 
the sign of the imaginary part of a is also changed as well. Then, by subtraction of these two 
Green functions, we get a delta function of v — a to find 



2ni 




PtW-Pl) 



(74) 



by omitting the pole parts. From this expression, eigenfunctions for E 
read 

Vc,kW) Y 7rsinh7rd/c 2 -p 2 \kTP 



k 2 /2 immediately 



(75) 
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Corresponding eigenfunctions in coordinate representation are also obtained by Fourier trans- 
formation as 



41w = 2 




9{±x)e lkr rF(l - ia,2,-2ikr), r 



\x\. 



(76) 



In the same way, eigenfunctions for bound states, given by Eq. ([7T|) , can be transformed to those 
of coordinate representation to yield 



Since a is not suitable for description of the Coulomb system, we may consider linear combi- 
nations of wave functions above for construction of eigenfunctions with even and odd parities. 
This completes the process of solving the Coulomb system in terms of path integral. 

To close this section let us add a comment on the difference in the number of bound states 
of the Coulomb system from that of H§. We have observed that, through the procedure of 
constructing (Hq — E) from the eigenfunctions of -ffs with the aid of uj — ip and its inverse, 
one of the eigenvectors of Hs has been eliminated from the series in Eq. ([53"]) . Therefore the 
dimension of the space spanned by the bound states of the Coulomb system is fewer by this 
missing one than that of Hs with E — — lu 2 /2. Since Hs is Hermitian, its eigenvectors form 
a complete set. Hence the bound states of He cannot be complete. This will become clear 
if we set a = 0, by which we return to the free particle, in Eq. (|67p . The sum over discrete 
eigenvalues is missing a term of n = 0, thereby failing to reproduce a delta function 5(j>b — p a )- 
Interestingly, the missing term is supplied from the term of continuous eigenvalues. The same 
can be seen more clearly in the integral of Eq.(|73p for the resolvent operator of a positive 
energy. For this case, if a ^ 0, poles from v/{e 27 ' v — 1) for bound states are complemented by 
the contribution from the pole of l/(v — a) for continuous eigenvalues to form a complete set. 
If we set a — after shifting v by —i/2(+i/2) for the first(second) term, the pole at v = a 
is absorbed into those of 1/(1 + e 27r!y ) by adding new poles at v = ±i/2. Cauchy's integral 
theorem then tells us that the sum over these poles at v — ±i(n + l/2)(n = 0, 1, 2, . . . ), with 
suitable regulators for each sum from upper and lower half-planes, is entirely equivalent to the 
integration of the continuous eigenfunctions over v along the real axis, therefore resulting in 
a delta function 5(pb — p a ) from both approach. It is exactly the pole at v = a in Eq. (|73|) 
that generates the term of continuous eigenvalues in Eq. (|67j) . Therefore, we are convinced the 
incompleteness of the set of eigenfunctions of bound states of the Coulomb system. It will be 
further interesting to notice that the eigenvalue, Eq = — oo, of the missing discrete eigenfunction 
is connected to that of E = k 2 /2 by letting k — > oo because, if a ^ 0, a — is possible only 
in this limit. Namely, the scattering state with an energy of positive infinity is complementary 
to the bound state with an energy with negative infinity though unphysical they are. A very 
intuitive understanding for this will be possible if we remember the behaviors of eigenfunctions 
in the coordinate representation; both limits, n — *■ so that uj n — > — oo for the bound state and 
k — > oo for the scattering state, yield a delta function S(x) although approaching in different 
directions. Finally, the absence of the bound state for n = in Eq. ([6"7f may be regarded as 
the consequence of the boundary condition imposed naturally by the formulation of the path 
integral for wave functions at x = 0. In view of the behavior of eigenfunctions at x = 0, we 
recognize that the path integral requires the Dirichlet condition at the origin. 

4 Conclusions 

On the basis of the Duru-Kleinert transformation with Fujikawa's gauge theoretical technique, 
we have formulated exact path integrals for the Coulomb potential in one space dimension. 
The mechanism of the exactness is explained by the pole structure of kernels in the momentum 
representation so that the time-sliced path integrals are evaluated simply by the use of Cauchy's 
integral theorem. Therefore the model considered here will be regarded as a new example of the 
path integral that can be performed exactly without making use of the Gaussian integration. 
It will be, however, difficult to observe such a structure if we employ the formal continuum 
formulation path integrals. A special care on the bound state with infinite binding energy has 



$J>(a;) = 2(-ir +1 io 3 n / 2 9(±x)e-^ r rF(l ~ n,2,2u n r), «, 



a 



\x\. 



(77) 
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been taken to clarify its absence. The disappearance of this unphysical state in the eigenfunc- 
tion expansion of the resolvent operator should be attributed to a proper prescription for the 
boundary condition at the origin required by the path integral formalism. 

The mechanism of the exact path integrals may have another interpretation; each relation 
between p a and pb read from the poles of the kernels in the momentum representation defines 
an action of the one parameter subgroup of SL(2, R) (setting r = it is assumed for the case of 
a negative energy) through a linear fractional transformation. This will bring us a group theo- 
retical approach for the understanding of the exactness. From this view point, generalizations 
of our technique to be useful for other systems will be interesting for future investigation. 
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